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(54) Title: CONCURRENT PROCESS FOR BLIND DECONVOLUTION OF DIGITAL SIGNALS 

^ (57) Abstract: This process has the objective to eliminate the problem of intersy mbol interference in digital signals, which is caused 
by the dispersive effect of any practical transmission channel. The solution for the problem of superposition of propagation rays 
O and its dynamic variation in the transmission channel of a digital system will be one of the hardest challenges for the technological 
improvement in this area for the next few years. In mobile communications, this undesirable improvement in this area for the next few 
J^J years. In mobile communications, this undesirable phenomena is characterized by parameters such as delay spread, angle spread and 
^ Doppler spread, which determine the level and dynamics of superposition between the system symbols, occurring as a consequence 
of the transmission of information through the channel. The devices that are likely to have their performance improved by the user of 
£5 the Concurrent Process for Blind Deconvolntion of Digital Signals includes but is not limited to: Spatial-Temporal Processing (used, 
for example, in Smart Antennas; Smart Sensors, etc...); Digital telecommunication systems in general (cellular telephony, digital 
^ television, digital radio, etc...); telemetry systems, remote sensing systems, geodesic localization/measurement systems (GPS, etc..), 
• "a storage systems, RADAR systems, 
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'"CQSCOKRENT PROCESS ICR BUND DBOIWCXimaVOFKGmSLSIGmLS". 

The concurrent process for blind deconvolution of digital signals described herein is 
intended to reduce and, possibly, eliminate the problem of intersymbol interference in digital 
signals, which is a result of the dispersive effect of any practical transmission channel. The 
5 process has the following working principle: The deconvolution is done by two filters, 
denominated herein V and W, whose two signals which results respectively from the 
filtering action of each filter jointly define the output signal y of the deconvolution process. 
The coefficients of filter V are adjusted by gradient methods, or equivalent, based on a cost 
function that measures a statistic dispersion of y. The coefficients of filter W are adjusted by 

10 gradient methods, or equivalent, based on a cost function that measures a distance fromy to 
the nearest digital alphabet symbol. The updating of the coefficients of filters W and V is 
such that the updating of the coefficients of W depends on the updating of the coefficients 
of V by means of a non-linear function. The results obtained demonstrate that the process 
described herein presents a performance comparatively superior to the processes usually 

1 5 adopted for the solution of this type of problem, and the range of possible applications and/or 
devices to be protected. Applications and/or devices subject to have their performance 
improved by the use of the Concurrent Process for Blind Deconvolution of Digital Signals 
includes but is not limited to: Spatial-Temporal processing (used, for example, in Smart 
Antennas or Smart Sensors, etc.), any digital telecommunication systems (cellular 

20 telephony, digital television, digital radio, etc.), telemetry systems; remote sensing systems; 
geodesic localization/measurement systems (GPS, etc); navigation aid systems; seismic 
survey systems by wave refraction/reflection, magnetic media storage systems, RADAR 
systems, SONAR systems. 

The Concurrent Process for Blind Deconvolution of Digital Signals, implemented by 

25 means of the Concurrent Equalizer, presents the capability to efficiently reach, under the 
same operational situation, a condition much closer to the zero-forcing condition, or 
equivalent, than the other equalizers used in the current scenario. This scenario includes the 
universally used - and renowned-CMA Equalizer. This allows the Concurrent Equalizer to 
be used to deconvolve channels with denser symbol constellation than the current 

30 technologic state permits. This way, the use of the Concurrent Equalizer for Blind 
Deconvolution of Digital Signals will allow for an increase in the transmission rate for those 
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digital systems that have their maximum speed limited by the operation of the equalizer (for 
example, see [21]). 

Definition of the Problem whose Solution is Obtained by means of the Concurrent 
Process for Blind Deconvolution of Digital Signals. 
5 The solution for the problem of superposition of propagation rays and the dynamic 

variation in the transmission channel of a digital system will be one of the hardest challenges 
for the technological advancement in this area for the next few years. In mobile 
communications, this undesirable phenomena is characterized by the parameters delay 
spread, angle spread and Doppler spread, which determine the level and dynamics of 

1CT superposition between the symbols of the system. Such symbol superposition occurs as a 
consequence of the transmission of information through the channel. 

Particularly, the superposition of symbols occurs due to the dispersive nature of the 
transmission channel of any practical digital system, which implies in the impulsive 
response e(n) of the channel to assume non-null values along many sample intervals [1]. A 

15 non-dispersive channel presents a response c(n) to the impulse S(n), characterized by a 
single S(n-n d ) impulse that occurs in n d ^ 0 sample instants after the initial instant. 

The dispersion of the channel may be better understood from the point of view of the 
convolution [2] between the sequence of symbols $(«) sent by the digital transmitter and the 
impulsive response c(«) of the channel. If the channel is non-dispersive (which is the ideal 

20 situation desired in practice, due to the high speed of transmission of the digital symbols of 
the current systems), then the symbol sequence «(») received by the digital receiver will be 
a equal to s(n), except for a gain constant G and a delay >0, i.e., u{n)=Gs(n-n d ).l£ 
the channel is dispersive, each sample of u(n) will be a weighted sum of the previous 
samples, which is defined by c(n) * 5 (n - « d ) . 

25 This way, the channel dispersion may be associated with the concept of 

reverberation, since the value of each sample is interfered by "echoes" from previously 
occurred samples. From this fact we have the concept of the quantity denominated Inter 
Symbol Interference or ISI, whose value measures how much the symbols of a digital 
transmission system overlap each other due to the channel dispersion [1][3]. 
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As the reverberation is a natural an inevitable phenomena to all and any practical 
transmission channel, the scientific community has circumvented the problem by means of 
dispersion compensation systems denominated Equalizers. Equalizers are usually placed 
between the digital receiver and the channel, so that the combined impulsive response h(n) 
of the Equalizer in series with the channel tends to S(n - n d ), and then avoiding the channel 
dispersion [1]. Figure 1 shows the simplified diagram of the dispersion compensation 
process accomplished by the Equalizer. 

m r 
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Transmission Channel 



Any media liable to 

transport Information by means of a group of 
digital symbols belonging to a digital alphabet 
For example, electro-magnetic, acoustic optical 
channels, etc... 



Process : The u(n) signal received from the channel is 
corrupted due to dispersion and/or echoes in the 
Equalizer input The Equalizer processes the u(n) 
signal and delivers the y(n) signal in its output, which 
presents contents of dispersion and/or echoes lower 
than the u(n) signal, due to the deconvolution process 
performed by the Equalizer. 



Figure 1: Baseband [1] representation of the removal process of the ISI by the Equalizer in 
a digital transmission system. 



In Figure 1, note that h(n) is the result of the convolution of the channel impulsive 
response c(n) with the Equalizer impulsive response /(«),i.e., h(ri) = c(ri)* /(«) [2], As 
the sequence of symbols >(«) received by the digital receiver is the result of the convolution 
of s(»)with h(n), and since the Equalizer actuates in the process in order that 
h(n)^S(p-n d ), the process of channel deconvolution - also called process of channel 
equalization [2][3] - results in y(n) -+s(n-n d ). That is, the deconvolution process 
performed by the Equalizer causes the y(n), contrary to the u(n) signal received from the 
channel, to be a reply of the s(ri) sequence originally transmitted, except for a delay of n d 
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samples. Every time the deconvolution process reaches a stage in which h(n) ->• S(n - n d ) it is 
said that the Equalizer has reached the ZF (zero forcing) condition (or situation), when then, 
ISI-»0 [1]. In many practical situations, due to several factors, such as noises and 
characteristics specific to the channel, the ZF condition is not fully accomplished, but only 
approximated by the Equalizer. One of these cases is that one in which the Equalizer operates 
in a situation of MMSE (Minimum Mean Squared Error), which is not exactly the ZF 
situation, but for practical and operational purposes of an Equalizer operating in this condition, 
it is a substitute situation. 

Independently of the transmission channel media (electro-magnetic, acoustic, etc.), 
due to the tendency of digital systems to transmit digital symbols at faster and faster speeds 
to meet market demands, the dispersion, even at low levels, is enough to deteriorate the 
system performance, because the duration interval of the symbols is short if compared to the 
duration interval of the dispersion. Another additional factor has recently appeared, 
compUcating even more the advancement of digital communications: The 
telecommunication market has presented a high growth rate of demand for more complexity 
(and, therefore, volume) of information to be transmitted. The Cellular Telephony System of 
3rd Generation based on the "software radio" concept (under development by several 
manufacturers [4]) is a typical example of such demand. The fact is that this demand growth 
rate is a lot higher than the operation speed growth rate (clock) of the digital circuits that 
compose the digital system, and that speed is limited by the status of technological advance 
in the micro-electronics and semiconductors field. Also, it is important to point out that all 
these limiting factors occur in a context with little availability to allocate spectral bands, 
especially in low frequencies. 

This way, in order to meet this higher demand for information volume, the solution 
found is to increase the 'Vocabulary" of digital symbols of a system (denominated 
constellation or system alphabet), so that the higher complexity of information may be offset 
by a better elaborated "vocabulary", capable to represent a maximum volume of information 
within the same time interval. That is, a denser constellation allows the transmission speed 
not to be reduced as a consequence of the higher complexity - i.e., entropy [1][3] - of 
information, since the larger number of symbols available circumvents the need to issue 
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many symbols in a row, a situation that would occur for a lesser dense constellation 

However, for the same number of "letters" available, a forcibly long 
"vocabulary" presents a larger number of "words" or symbols that are similar to each other 
than a simpler "vocabulary". This increased similarity between the constellation symbols of 
a digital system that uses a denser constellation turns into a serious problem for the system 
performance when the channel is dispersive. This is due to the fact that the minimum ISI 
makes the symbols unintelligible from the digital point of view, due to the increased intrinsic 
similarity between them. In technical terms, the distance that differs the symbols from each 
other becomes shorter, which hampers the identification of the symbol effectively received. 
This way, the performance of a channel Equalizer of a digital system becomes the critical 
factor in the performance of the system as a whole. 

It is important to point out that the emerging cellular wideband 
communications systems, such as the 3rd Generation systems, involve the use of Smart 
Antennas in radio-bases and, possibly, in receiver terminals. These Smart Antennas are 
controlled by digital processing algorithms, denominated "Spatial-Temporal" algorithms, 
which are essentially equivalent to the equalization algorithms referred to and/or described 
herein. 

An equalizer is said to be blind when it performs the blind deconvolution of the 
channel, that is, when the equalizer totally discards the need of the receiver to recognize the 
s(n) sequence originally transmitted. On the opposite, non-blind equalizers need tixQs(n) 
sequence to be recognized in the receiver as a reference - at best for a short interval of time 
- so that the ZF, MMSE or other equivalent condition may be accomplished. The advantage 
of a blind equalizer over a non-blind one is, therefore, obvious from the point of view of the 
reduction of spectral band required. 

The Concurrent Process for Blind Deconvolution of Digital Signals, objective of this 
patent claim, is implemented by means of a blind equalizer system whose architecture is 
characterized for being concurrent. Therefore, hereinafter, we will refer to the equalizer that 
implements the Concurrent Process for Blind Deconvolution of Digital Signals as 
"Concurrent Equalizer", The new and inedited architecture of the Concurrent Equalizer - 
and the operation mode intrinsic to this architecture - permits the efficient deconvolution of 
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digital signals with much denser constellations than permitted by the types of blind 
equalizers usually employed in current digital systems. Furthermore, the Concurrent 
Equalizer permits the dynamic temporal variation in the channel, resulting, for example, 
from the relative displacement between digital transmitter and receiver, to be significantly 
5 higher than the maximum displacement speed usually permitted by other equalizers. 
The Concurrent Equalizer. 

Figure 2 shows the diagram for the Concurrent Equalizer. From its initialization, the 
vector with weights V = [V 0 V l A V LA J is updated with some algorithm based on 
__ Gradient methods, such as, for example, Stochastic Gradient [3], RLS (Recursive Least 

10 Square) [3], or another equivalent, intending to minimize the cost function J Dl of another cost 
function that measures dispersion. Likewise, from its initialization, the vector with weights 
W = [W 0 A W M J. is updated with some algorithm based on Gradient Methods, 
intending to minimize the cost function Jq, or another that measures an output distance y 
from the Equalizer to the nearest symbol in alphabet A, which is given by Q{y}. The 

15 operator Q{} represents the quantization operation with reference to the symbols of 
alphabet A . 
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QM = argmin|y-5 k | ,s a6 a 
A =0,1, ...,M-1 




Stochastic Gradient 

Figure 2: Schematic diagram of the Concurrent Equalizer, which implements the Concurrent 
Process for Blind Deconvolution of Digital Signals. Each z" 1 block introduces the delay of a 
sample. 

5 In Figure 2 y = e{|A| 4 }/e{| A| 2 J is the statistic dispersion constant for a dispersion order of 
P=2. E{} is the operator that returns the statistic average of the argument. The {f operator 
results in the transposition of the vector/matrix argument and the j • | operator returns the 
Euclidean rule of the argument. The | • | 2 operator was used as a distance metric, since the 
symbols s of alphabet A are usually complex numbers (when, for instance, the signalling is 

10 M-QAM, M-PSK, [1] etc.). Therefore, the computational cost of the |.| 2 operation is 
reduced with | • | 2 = s • s' , where the {■}* operator represents the conjugation of the argument. 
However, any operator expressing the length - in the distance context - of the argument may 
be used. The Gradient method shown in Figure 2 is the Stochastic Gradient, but any 
equivalent method is liable to be employed. 

15 Note that J D is a statistic dispersion of order P [5], which inherently uses higher 

order statistics from the equalizer output set y [6]. Although the Concurrent Equalizer can 
operate for any internal value of P, for the sake of simplification of the following 
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description, we will adopt P = 2 . 

Note that Jq measures the square of the Euclidian rule between the equalizer output 
and the nearest alphabet symbol. Also, by the definition of the cost functions, note that if 
Jq is minimized, Jd will simultaneously be minimized too. However, the reciprocal is not 
5 true. 

In order to the nnnimization process of Jq does not destroy the stability of the 
delicate dynamics of the minimization process of Jd, for both share the same common 
variable;/, a link is established between them, as demonstrated in Figure 2. This link permits 
the process of J D to inform how balanced it is to the process of Jq. The balance of the J D 

10 process must be understood not only associated with the stability concept, but also associated 
with the degree of the effect resulting from that stability, that is, the degree in which the 
minimization process of Jd simultaneously minimizes Jq, since the equalizer was developed 
to operated concurrently. This criteria is critical, since the minimization of Jd, for definition, 
not necessarily minimizes Jq. This link consists of a non-linear operator that inhibits the Jq 

15 process when the minimization process of Jd does not minimize Jq simultaneously. Thus, 
although they share the common output variable y - which makes the processes concurrent 
due to the resulting architecture - the Jq process can only interfere cooperatively with the 
minimization process of Jd. 

Table 1 describes the general procedure adopted when the Concurrent Equalizer is 



20 applied to the equalization of a channel under T/2 fraction sampling [1]. 



Stage 


Procedure 




Initialize vectors W and V: 

W = 0+j0 

and 

v (0+j0, *=OMW, k*€ 
n -\l+j0,k = ? 

where, j = >PT , L is the dimension of vectors V and (equalizer size) and £ , 
0<£<I-1, is the index of the only weight of the V weight vectors to be 
initialized with the value 1+jO. 




Initialize the index counter for the received samples by fractional sampling from 
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channel: z' = l 




Initialize the channel regressor index counter(or instanl/iteration index counter): 
n=0 




Obtain the rP channel regressor r{n) [6][1 ][!]: 
n(n)=u(L-l-k+t), * = 0,1,A,L-1, 

T 

where, v is the sequence of samples received by fractional sampling ~ from the 

channel, with i = 1,3, A ,N a - 1 , varying as » = 0,1, A ,N r - 1 , and is the total 
number of samples to be received by fractional sampling from the channel, 

jVr=|^— is the total number of regressors to be obtained from the 

channel and T is the interval between the symbols generated in the transmitter. [-J 
is the operator that results in the closest whole number minor than the argument. 




Obtains the equalizer output at the instant n: 
y(n) = W T (n)-r(?7)+V r (nyrM 




Update V weight vector: 

V(n+l) = V(n)+tj v .y(n)( r -\y(nf)-r\n), 

where, tj v is the adaptation step of vector V, 0 < % « 1 .0 . 




Update W weight vector: 

E(»+D=E(")+^[i-^(»)][Qk»)}-X#(») . 

where, rj w is the adaptation step of vector W, 0 < i] w « 1.0 . 

fo, Q{y(")}=Q{y(»)} 
A>(»H r i r -i 

controls the non-linear link of Figure 2, where 
y(«) = r r (n + 1) • r(n) + W T (n) ■ r(n) 




Increment index counters: 

i=i+2 

«=«+l 




Test loop end: 

If I +/>#„=> END 
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[ | Otherwise, repeat steps 4 to 9. ' I 

Table 1: General summary of operations performed by the Concurrent 
Equalizer to implement the Concurrent Process for Blind Deconvolution of Digital Signals 
when applied to the channel equalization, for p=2. Note that, in a continuous 
communication, the value of N a may be considered infinite. 
> From Table 1, note that factor [l-D B ] generates the following behavior during the 

operation of the Concurrent Equalizer: W can only be updated if, when updating £ the 
output variable y remains in the same quantization level Q{y} . A quantization in V that does 
- not change the quantization level Q implicitly means a probably correct updating, since this 
is the expected condition after the equalizer convergence. Therefore, factor [l-D s ] controls 
' the updating of W pursuantly to the concurrent operation of the system: The minimization 
process of J Q is only activated when the niinimization process of J D simultaneously 
mimmizes J Q . As the E and if filters jointly influence the output* factor [l-D Q ] avoids 
any non-cooperative interference of the minimization process of J Q that might unbalance the 
minimization process of J D . 

Comparative Tests under Static Channel - Results Obtained 
With the objective of testing the Concurrent Equalizer in situations as real as 
possible, we will use in this description models of channels available in the database of Rice 
University, in Houston, Texas, USA This database is known as Signal Processing 
Information Base (SPB) [8] and can be accessed at http://spib.rice.edu/. The SPIB has been 
widely used as a reference in many recent works that intend to avoid an exclusively 
academic focus. Additional information, including studies and articles based on SPIB can be 
found in the site of the Blind Equalization Research Group (BERG) [9] at 
http://backhoe.ee.cornell.edu^ERG/, which belongs to the Cornell University, in Ithaca, 
New York, USA The operation of the Concurrent Equalizer under dynamic channel is 
described in Section 4. 

The microwave channel models to be used may be accessed at 
http://spib.rice.edu/spib/microwave.html. These models constitute the impulse response of 
several real microwave channels obtained from field measurements. Although it is herein 
presented the complete characterization of each channel in time domain and also in 
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frequency domain, additional information may be obtained in [10], [1 1] and [12]. 

The impulse response of the SPIB microwave channels is obtained under a 
high sampling rate, at the magnitude of tens of megabauds per second, resulting in impulsive 
responses with hundreds of samples. This allows each researcher to decimate the resulting 
5 sequence by a factor appropriate to its particular interest, without losing significant 
information. The great majority of the works rarely uses an impulsive response of channel 
with more than some tens of samples. In this description, in order to maintain a comparison 
term with recent works in the field [6][7][13], the SPEB microwave channels will be 
decimated to 16 samples. Table 2 identifies the SPIB channels used and establishes the 
10 respective designation to be adopted hereinafter. For example, references to channel M4 in 
this description must be understood as: the channel whose impulse response results in the 
decimation to 16 samples from the sample sequence of the file chan4.mat from the SPIB 
database, concerning microwave channels. 



SPIB File: 
(microwaves) 


Sampling rate l/T of 
transmitter in Mbauds/s 
(according to 

http://spib.rice.edu/spib/toic 
rowavehtml): 


Number of samples in the 
SPIB file, resulting from the 
fractional sampling T/2 of 
impulsive response received 


Designation, in this 1 
description, of the respective 
impulsive response generated 
by decimation of the SPIB file 
to 16 samples: 


chan4.mat 


30 


300 


M4 


chan9.mat 


30 


300 


M9 


chanl4.mat 


30 


600 


M14 



Table 2: Microwave channel models in the SPIB database. 



15 In order not to change the location of the zeros of the FIR transfer function of 

a SPIB channel, the decimation is applied to the frequency domain [2] [7]. Specifically, the 
number of samples in a SPIB file is extended to the next whole power of 2 plus the 
addition of null samples in sequence. Then, the FFT (FFT - Fast Fourier Transform) is 
applied on the extended sequence and the resulting sequence in the frequency domain is 

20 decimated to 16 samples. Next, the IFFT (IFFT - Inverse Fast Fourier Transform) [2] is 
applied over the decimated sequence in the frequency domain, resulting in a sequence of 
16 samples in the time domain, corresponding to the decimated impulsive response. If the 
decimation in the time domain were used, aliasing effects would be generated [2], 
changing, particularly, the position of reflected roods and roots close to the unitary circle 

25 in the frequency domain z, and then distorting the characteristics of the SPIB channels. 
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Another channel model used in this description is the representative of cable TV 
channels, which may be accessed at http://spib.rice.edu/spib/cable.html. The sampling rate 
used for these SPIB channels is the usually employed in practical cable TV decoders. This 
way, in order to maintain a comparison term with the works in the field, these channels will 
not be decimated. Table 3 identifies the cable TV channel SPIB to be used and establishes 
the respective designation to be adopted hereinafter. For example, references to channel CI 
in this description must be understood as: the channel whose impulse response is given by 
the sequence of samples of the chanl.mat file of the SPIB database concerning cable TV 
channels. 



SPIB File: 
(cable TV) 


Number of sample sin the SPIB file, 
resulting from the fractional sampling J/2 
of the impulsive response received from the 
channel: 

(according to 

http://spib.rice. edu/spib/cable.html.) 


Designation in this 
description of the 
respective impulsive 
response associated: 


Chanl.mat 


128 


CI 



Table 3: Model of cable TV channel from the SPIB database. 



In this description, the channels M4, M9, M14 and CI are characterized by: 

1- Module [c(n)| and phase Zc(n) (in degrees) of the impulsive response c(n) resulting 
from the fractional sampling 7/2 of the impulse response c(r) of the channel. 

2- The location in plane z of the zeros of the transfer function C (z) of the channel. 

3- The place in plane z of the zeros of the transfer functions C EVEN (z) and C 0DD (z) of the 
sub-channels (fractional sampling T/2 implies odd and even sub-channels [6][7]). 
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Figure 3: Characterization of channel M4. and 
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Figure 4: Characterization of channel M9. 
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Fig. 5c: Zeros of C(z) . Fig. 5d: Zeros of c^oo - A and c 000 ^) ~ * 

Figure 5: Characterization of channel Ml 4. 
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Figure 6: Characterization of channel CI. 
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The Concurrent Equalizer performance evaluation will be done in reference with an 
equalizer type CMA (Constant Modulus Algorithm) [1][7]. The CMA FIR transversal 
Equalizer is perhaps the most widely used and tested in the implementation of the current 
blind equalizers [14], because, in practical systems, its performance has been reported to be 
5 superior to other blind equalizers currently in use [6] [7]. The performance criteria adopted 
are based on the following simulation results: 

1- Graph of the constellation T at the output y of the CMA Equalizer after convergence, 
having the constellation of the source alphabet A as reference, with A represented in the 
graph by the symbols "+ ". 
10 2- Graph of the *F constellation in the output y of the Concurrent Equalizer after 
convergence, having the constellation of the source alphabet A as reference, with A 
represented in the graph by the symbols "+ ". 

3- Comparative graph of ISI(/) curves for the two equalizers, i - 0,2,A ,N a -1 , where N a is 
the total number of samples received by fractional sampling T/2 from the channel, and 

15 ISI(0 is the "reverberation" measure in the equalized signal y in instant /', obtained from 
ISI= [Z[ff*| - n^ff t |y ma^l calculated for this instant. H k represents the ft* component of 
vector H, which contains the coefficients of the impulse response h(n), combined with the 
equalizer and the channel. 

4- Comparative graph of MSE(7) curves for the two equalizers, / = 0,2,A ,N a -l, where, 
20 MSE(z) is the mean square error between the equalized output y and the original sequence 

transmitted s, taking into consideration for the average the interval of samples that includes 
the instant i and the L+L,-l samples previous to it. L c is the dimension of channel 
dispersion defined by the number of non-null samples in the impulse response c(n) of the 
channel. Additionally, this graph demonstrates the horizontal line NT, which defines the 
25 maximum permitted MSB level, so that the transference to the renowned DD Equalizer 
(Direct Decision) [1][6][7] be done. The DD Equalizer acts after the convergence of the 
CMA Equalizer, intending to reduce its residual MSE (steady state MSE). However, the DD 
Equalizer can only converge if the steady state MSE of the CMA Equalizer has reached a 
value inferior to the indicated by the horizontal line NT [7]. 
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The sequence representing the information source s to be transmitted through the 
channel has average zero, unit variance, uniform distribution, and presents statistic 
independence between samples. 

The noise added to the equalizer input is Gaussian, with zero average, spectrally 
5 white and with variance defined by SNR (Signal To Noise Ratio) [ 1 ] . 

Some authors suggest that the dimension of equalizer L be done larger than the 
dimension of the channel dispersion L c , so that the extra degrees of freedom so obtained 
allow the rninimization of the equalizer noise gain [15]. However, in order to minimize the 
hardware cost for implementation purposes, we will adopt L = L e in the following 
10 simulations. 

The initialization of the vector that defines the filter of the CMA Equalizer and the 
vector V of the Concurrent Equalizer follows the strategy denominated single spike [7], which 
consists of imtializing the vectors with l+_/0 in position 0<| <Z-1, rnamtaining the 
elements in other positions null. In [7] it is demonstrated that the optimal value of £ depends 

15 on the position of the sample with maximum absolute value in the channel impulse response 
\c(n)\ sequence, denominated channel cursor. However, as c(«) is not known at first, the 
universally accepted practice is to make £=Z/2. If the equalizer does not converge with 
£ = L/2, it is normal that the receiver system tries to make a few attempts around £ = L/2 
until the equalizer converges, that is, g=l/2±A, A = l,2,A [16]. We will adopt such 

20 heuristic in case of initial convergence failure. The cases for which it was necessary to apply it 
can be identified by £ * L/2 . 

The adaptation step of the CMA Equalizer filter is n, and the adaptation steps of the 
Concurrent Equalizer filters are tj v and tj w . These parameters are experimentally adjusted in 
order to meet the tradeoff between fast convergence and low steady state MSE. 

25 With unitary variance in the information source s, for 16-QAM and 64-QAM, tj 

presents better results in the neighborhood of value lxlO" 3 , and for 256-QAM, in the 
neighborhood of lxlO" 4 . These results agree with the experiments involving 16-QAM and 
256-QAM modulation in [7]. 
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In the following graphs, the foEowing events are considered indication of good performance: 

1- The degree in which the points of the output constellation y of the equalizer (TorW) 
concentrate around the + symbols of alphabet A. Under ideal performance, this graph 
would show all of the output points y projected in the intersection of the lines that define 

5 the + symbols and the alphabet A . 

2- The rate at which the curve ISI(/) decreases with i and how low and invariable is the 
value of ISI(i) after convergence (steady state ISI). 

3- The rate at which the curve MSE(i) decreases with i and how low and invariable is the 
— value of MSE(/) after convergence (steady state MSB). 

10 Below are the results of the simulations which assess the comparative performance 

between the CMA Equalizer and the Concurrent Equalizer. 
Modulation 16-QAM SNR=35dB - Channel M14 




Fig.7a: T (CMA) Fig. 7b: W (CONC) 



WO 02/071607 



19 



PCT/BR02/00030 



Fig. 7c: ism 



Fig. 7d: MSE(/) 



Figure 7: Comparative results between the CMA Equalizer and the Concurrent 
Equalizer (CONC) for channel M14, modulation 16-QAM (^ = 1.32), SNR=35dB. 
Parameters: N a =100000, 77 = 1 xl0~\ /^lxlO -3 , rj w = lxl0'\ L = L C =16 and £ = 6. 
Modulation 64-QAM SNR=35dB - Channel M9 




Fig. 8a: T(CMA) 



Fig. 8b: ¥(CONC) 
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Fig. 8c: ISI(i) 



Fig. 8d: MSE(i) 



Figure 8: Comparative results between the CMA Equalizer and the Concurrent 
Equalizer (CONC) for channel M9, modulation 64-QAM =1.380953), SNR=35dB. 
Parameters: N a = 100000 , 77 = lx 10" 3 , tj v = lxlO" 3 , tj* = lxlO" 2 , L = U = 16 and £ = 2 . 
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Modulation 256-QAM SNR=35dB - Channel CI 
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Fig. 9a: T(CMA) 



Fig. 9b: ^(CONC) 
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Fig. 9c: ISI(;) 
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Figure 9: Comparative results between the CMA Equalizer and the Concurrent 
Equalizer (CONC) for channel CI, modulation 256-QAM (y = 1.395295), SNR=35dB. 
Parameters: N a =100000, 77=lxl0" 4 , tj v =1x10" 4 , 77 w = 1x10~ 3 , L = L e =128 and £ = 63. 
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Figure 10: Comparative results between the CMA Equalizer and the Concurrent 
Equalizer (CONC) for channel M4, modulation 16-QAM (/ = 1.32), SNR=35dB. 
Parameters: 2V fl =30000, ^lxlO" 3 , ? 7v = lxlO~ 3 1 77^ =lxl0" 2 , L = L C =16 and £=8. 
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Modulation 64-QAM Sm=35dB - Channel M4 
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Fig. 11a: T(CMA) 




Fig. lib: T(CONC) 




Fig. 11c: ISI(0 
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Figure 11: Comparative results between the CMA Equalizer and the Concurrent 
Equalizer (CONC) for channel M4, modulation 64-QAM 0 = 1.380953), £NR=35dB. 
Parameters: ^=30000, J7 = lxl0" 3 , t7 v = lxl<r 3 , J7 w =lxl0~ 2 , L=L C = \6 and £=6. 

Note that, in none of the first three examples shown (Figures 7, 8 and 9), the 
CMA Equalizer reaches a MSE low enough so that the transference to the DD Equalizer is 
successfully accomplished. These are typical cases in which the popular and renowned 
CMA-DD Equalizer would fail in meeting its objective. For the three examples, the 
Concurrent Equalizer not only has surpassed the MSE level for transference, but also 
reached a MSE level much inferior and less variable than the CMA Equalizer. The same 
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observation is valid for the comparative ISI between the equalizers. 

The Concurrent Equalizer outperforms the CMA Equalizer, even in those 
situations in which the CMA reaches the transference MSE, such as demonstrated, for 
example, in Figure 10 for modulation 16-CAM. Note also that for the same channel and 
5 modulation 64-QAM, the CMA Equalizer reaches only a marginal steady state MSE, as 
shown in Figure 11. 

Note that in the examples for which the channel rotates the received 
constellation, the Concurrent Equalizer is capable to correct this undesirable effect, while the 
CMA Equalizer does not have this capability. This is a particularity of the Concurrent 

10 Equalizer, which allows the simplification of the synchronism system. It is important to 
point out that, upon performing joint equalization with carrier phase retrieval and/or symbol 
sampling instant retrieval, the Concurrent Equalizer presents an inertia and a response time 
shorter than those of the CMA Equalizer followed by a phase corrector (de-spinner). This is 
a decisive characteristic for the good performance of the Concurrent Equalizer in time 

15 varying channels. 

It is important to reaffirm that the Concurrent Equalizer surpasses the 
performance of the CMA Equalizer, not only regarding the MSE, but also the ISI. In all 
cases, the Concurrent Equalizer quickly reached a much lower and more stable ISI level than 
CMA Equalizer. This behavior of the Concurrent Equalizer demonstrates that its best 

20 performance compared to the CMA Equalizer does not depend on the phase rotation 
imposed by the channel to the constellation received, since the ISI is defined from the 
module of the impulse response h(ri). If, in the examples above, a system to correct the 
phase rotation on the received constellation were previously added to both equalizers, the 
performance of CMA Equalizer would not be changed regarding the ISI due to the phase 

25 invariance of this equalizer [7], while, on the other hand, the performance of the Concurrent 
Equalizer would tend to improve. 

Comparative Tests under Dynamic Chatmel- Results Obtained 

In the following simulation set, the performance of the Concurrent Equalizer is 
evaluated against the CMA Equalizer, both operating under temporal variation in the 

3 0 channel situation which characterizes the operation under dynamic channel. 

The objective is to determine the sensibility of both equalizers to the speed of 



WO 02/071607 



24 



PCT/BR02/00030 



variation in the channel. To do so, the following analysis is based on the channels suggested 
by M. Gosh [17], whose work referred to herein studies the digital television system ATSC 
8-VSB proposed by the Advanced Television System Committee [18] for the equalization 
performance. 

5 The ATSC 8-VSB digital system was proposed in the USA in 1993, with the 

objective to substitute the veteran NTSC analog system for television broadcast. Using the 
same channel with the 6MHz bandwidth of the NTSC system, the ATSC 8-VSB was 
conceived to present a superior performance, even under high levels of interference, noise 
and multi-path effects [19]. 

10 Since this description considers all signals as baseband signals[l], for the sake 

of this simulation, it is enough to use as information source s the 8-VSB alphabet generator 
A={-1.53 -1.09 -0.66 -0.22, 0.22, 0.66, 1.09, 1.53} to represent the eight possible signal 
amplitudes (8-PAM) baseband generated by a 8-VSB transmitter. Since the Channel Encoder 
of a ATSC 8-VSB transmitter "randomizes" with uniform distribution the baseband signal 

15 for transmission efficiency purposes [20], the source s presents a uniform probabilistic 
distribution with statistic independence between samples. Note also, from the definition of 
alphabet A , that the variance of s is unitary. The unitary variance is assumed for s as partial 
representation of the CAG system effect present in any digital system. 

The dynamic channel model adopted in this study is the one suggested in [17] 

20 and is specified in Table 4. 



Components of the 
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Delay [samples] 


Principal Ray: 
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Table 4: Specification of the dynamic channel mode adopted. The column 'Delay 
[samples]" is obtained by dividing the values in column "Delay [^s]" by the fractional 
sampling interval T/2 of the receiver, where \/T =F, =10. 1 '6MHz is the sampling rate of 
the symbols in the ATSC 8-VSB transmitter [19][18], considering as a sample of index 0 
5 the one with lower delay. The integer variable n represents the index of the n* channel 
regressor received, or the index of the n tt symbol sent by the transmitter and, therefore, 
presents temporal spacing T associated to the interval between the indexed samples. The 
real constant / represents the speed of temporal variation for each type of dynamic channel 
_ to be considered. For example, 0.05Hz < f < 0.5Hz models dynamic channels resulting 
10 from the reflection of signals by a tall building whose structure is swaying in its natural 
frequency of mechanical oscillation (due to the wind, etc.), and 2Hz<f <5Hz models 
channels with rapid and periodic fading (flutter) generated by the reflection of the signal by 
aircrafts flying in low altitude [17]. 

Therefore, from Table 4, the impulsive response obtained by fractional sampling 
15 T/2 of the dynamic channel in question has its representation in the z domain given by 

C(z) = 0.1+1.0z- 39 +0.1z" 42 +0.316 sin^/j^jz^+O^'+O.^z- 426 . (1) 

The ATSC 8-VSB system includes a baseband filter (Nyquist filter) of the 
raised-cosine type, with a roll-off of 11.5% [18][17]. As this filter affects the baseband 
signal, the channel model to be considered must include the effect of this filter in the 
performance evaluation of the equalizer. The raised-cosine filter of the ATSC 8-VSB system 
20 is an analog filter that presents a transfer function H(f) , according to Figure 12. 
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Figure 12: Graphic representation of the transference function H(f) of the 
raised-cosine filter of the ATSC 8-VSB system. 

The discrete impulsive response obtained by fractional sampling T/2 of the 
continuous impulsive response of the analog filter defined by the graph of Figure 12 is given 
by Equation (2) and is graphically represented in Figure 13. 



-(•HI 



(2) 



In the Equation (2) N=64 is the number of samples considered in this study for the 
impulsive response h n (/') , and a = 0, 1 15 is the parameter that defines the 1 1 .5% roll-off of 
the filter in the frequency domain. The integer variable i represents the index of the r* 
10 sample received from the channel and, therefore, presents temporal spacing T/2 i 
with the interval between the indexed samples. 
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Figure 13: Graph of h n (i) for N=64 a?ida = 0.U5. 

In the following simulation set, the effect of the raised-cosine filter is included by 
— effecting the convolution of the impulsive response of the channel defined by (1) with the 
5 impulsive response of the raised-cosine filter given by (2). Therefore, the dimension L e of 
the global channel dispersion ? including the effect of the raised-cosine filter ? is given by 
the sum of the dimension of channel dispersion with the dimension N of the filter dispersion, 
that is, h =427+64 = 491. 

The performance criteria adopted are based on the following simulation results: 

10 1- Graph of the CMA Equalizer output yauQi), where n is the index of the n th channel 
regressor received, having the source alphabet constellation A as reference. A is 
represented in the graph by the set of dashed lines. 
2- Graph of the Concurrent Equalizer output ;w(h), having the source alphabet 
constellation A as reference. A is represented in the graph by the set of dashed lines. 

15 3- Graph of the MSEau{n) curve of the CMA Equalizer, with MSEauin) being the mean 
square error between the yaa(ji) output and S{ya«(«)}, Q{} is the operator that 
represents the 8-VSB quantizer. The averaging interval for the MSEau(n) includes the 
instant n and the previous L c -1 samples. Additionally, this graph shows the horizontal 
line NT that defines the maximum permitted MSE level, so that the transference to the 

20 DD Equalizer be possible. 

4- Graph of the MSEcanc(n) curve of the Concurrent Equalizer, with MSEcaNc(ri) being the 
mean square error between the yooNcin) output and (2{)W(>j)}. The averaging interval 
for the MSEccnc(n) includes the instant n and the previous L e -1 samples. Additionally, 
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this graph shows the horizontal line NT that defines the maximum permitted MSE level, 

so that the transference to the DD Equalizer be possible. 

A SNR= lOOdB was adopted, since the objective is to measure the sensibility of the 
equalizer to the channel variations, and not the sensibility to the noise. As the channel varies, 
it was considered reasonable to adopt a dimension of the equalizer L that exceeds in 50% the 
dimension L c =491 of the global channel dispersion, that is, it was adopted Z = 738. The 
initialization of vector V, both for the CMA Equalizer and the Concurrent Equalizer, 
follows the usual technique single spike, which, in this simulation, consists of initializing V 
with the value 1.0 in position £ = 65 , maintaining null the elements in the other positions. 
The adaptation steps were fixed in 77 = 2 x 1 0 -4 for the CMA Equalizer and tj v = 2 x 1 0" 4 and 
rj„ = 2xl0~ 3 for the Concurrent Equalizer. 

Below are the simulations of the comparative performance between the CMA 
Equalizer and the Concurrent Equalizer, both operating under dynamic channel for variation 
speeds f = 0Hz (static channel), f = 10Hz, f=20Hz, f = 50Hz and f = 100Hz . 
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Fig. 14c: MSEo«(w) Fig. 14d: MSE«ac(») 

Figure 14: Comparative results between the CMA Equalizer and the Concurrent Equalizer 
(CONC) for the channel specified in Table 4 with f~0Hz, including the effect of 
raised-cosine filter. 
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Figure 15: Comparative results between the CMA Equalizer and the Concurrent Equalizer 
(CONC) for the channel specified in Table 4 with f = \QHz, including the effect of 
raised-cosine filter. 
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Figure 16: Comparative results between the CMA Equalizer and the Concurrent 
(CONC) for the channel specified in Table 4 with / = IQHz , including the 
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Figure 17: Comparative results between the CMA Equalizer and the Concurrent Equalizer 
(CONC) for the channel specified in Table 4 with f=50Hz, including the effect of 
raised-cosine filter. 
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Figure 18: Comparative results between the CMA Equalizer and the Concurrent Equalizer 
(CONC) for the channel specified in Table 4 with / -\00Hz, including the effect of 
raised-cosine filter. 



The experimental results under dynamic channel operation demonstrate that the 
Concurrent Equalizer presents a performance superior to the CMA Equalizer in the task to 
deconvolve channels that present rapid temporal variations. 
5 Note that, opposite to the Concurrent Equalizer, for / = 50Hz and / = lOOife 

the CMA Equalizer cannot reach and maintain the MSE lever required for the transference to 
the DD Equalizer. Note also that, for all cases, the MSE regime level reached by the 
Concurrent Equalizer is much lower and more invariable when compared to the level 
accomplished by the CMA Equalizer. 
10 These results suggest that the Concurrent Equalizer will be able to largely 

contribute for the performance improvement of mobile digital systems. 
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CLAIM 

1 "CONCURRENT PROCESS FOR BLIND DECONVOLUTTON OF 
DIGITAL SIGNALS", characterized in accordance with the descriptive report, in which it is 
supported by schematic representations of the process, where this Invention Patent claims the 
rights on (1) any application and/or device that uses a system for deconvolution of digital 
signals and that, to do so, is based on the following working principle: the deconvolution is 
done by two filters, denominated herein V and W, whose two signals which results 
respectively from the filtering action of each filter jointly define the output signal y of the 
deconvolutor system; the coefficients of filter V are adjusted through gradient methods based 
on a cost function that measures a statistic dispersion of y; the coefficients of filter W are 
adjusted through gradient methods based on a cost function, denominated herein J Q , which 
measures the distance from y to the nearest digital alphabet symbol, this symbol is 
denominated herein <2{y}, where the operator £){•} represents quantization; the updating of 
the coefficients of filters W and V is such that the updating of the coefficients of W depends 
on the updating of the coefficients of V by means of a non-linear function; (2) any application 
and/or device that uses a system for the deconvolution of digital signals and that, to do so, 
uses the working principle described in (1) for the operation of Decision Feedback 
Equalizers - DFEs; (3) any application and/or device that uses a system for the deconvolution 
of digital signals, whose operation is based in the working principle described in (1) and/or 
(2), in which Q{y}, required to the computation of Jq, is substituted with a temporal 
periodicity T during a time interval r , v <, T , by the value obtained from a framing sequence 
d sent by the digital signal transmitter; the training sequence d represents the sequence of 
symbols originally sent by the transmitter; the training sequence d is available to the 
deconvolutor system during the time interval T; the training sequence d is used as a 
reference for the deconvolutor system; (4) any application and/or device that uses a system 
for the deconvolution of digital signals, whose operation is based on the working principle 
described in (1) and/or (2) and/or (3), independently of the nature of the transmission channel 
physical medium (electromagnetic, acoustic, optical media, etc.); (5) any application and/or 
device that uses a system for the deconvolution of digital signals, whose operation is based on 
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the working principle described in (1) and/or (2) and/or (3 ) applied to Spatial-Temporal 
Processing (used, for example, in Smart Antennas or Smart Sensors), independently of the 
nature of the transmission channel physical medium; (6) any simulator, implemented in 
software and/or hardware, that uses a system for the deconvolution of digital signals, whose 

5 operation is based on the working principle described in (1) and/or (2) and/or (3); (7) the 

range of possible applications and/or devices to be protected, applications and/or devices 
subject to have their performance improved by the use of the Concurrent Process for Blind 
Deconvolution of Digital Signals includes but is not limited to: Spatial-Temporal processing 
_ (used, for example, in Smart Antennas or Smart Sensors, etc.), any digital 

.0 telecommunication systems (cellular telephony, digital television, digital radio, etc.), 

telemetry systems; remote sensing systems; geodesic localization/measurement systems (GPS, 
etc); navigation aid systems; seismic survey systems by wave refraction/reflection, magnetic 
media storage systems, RADAR systems, SONAR systems. 
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